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1  Introduction 


This  report  summarizes  the  results  obtained  in  a  three-year  investigation  into  the  application  of 
wavelet  methods  in  modulation  techniques  for  multi-user  spread  spectrum  communication  systems. 
The  investigation  was  undertaken  between  November  1992  and  January  1996. 

2  Objectives 

The  objectives  of  this  project  remained  essentially  as  outlined  in  the  project  proposal  and  pre- 
vious  status  reports: 

•  Task  1:  Development  of  and  evaluation  of  wavelet  symbols  for  multiple  access  spread  spectrum 
communications. 

•  Task  2:  Evaluation  of  communication  system  performance. 

•  Task  3:  Investigation  of  possible  countermeasures  to  wavelet-based  spread  spectrum. 

In  addition  to  addressing  these  objectives,  some  basic  research  on  wavelet  analysis  was  performed 
in  connection  with  the  project  which  led  to  a  class  of  multidimensional  generalizations  of  the  wavelet 
transform. 

3  Status  of  Effort 

The  project  was  completed  in  January,  1996  with  the  benefit  of  a  two-month  extension  of  the 
original  three-year  time  line. 

Task  1  led  to  development  and  analysis  of  a  class  of  bandlimited  orthonormal  wavelet  sets  having 
properties  weH  suited  for  use  with  the  wavelet-based  scale-division  multiple  access  communication 
scheme  around  which  much  of  this  project  was  based. 

Crosstalk  effects  in  scale-division  multiple  access  under  asynchronous  and  distributed  operation 
were  studied  under  Task  2.  Other  performance  questions  originally  planned  to  be  addressed  under 
this  Task  were  de-emphasized  in  later  stages  of  the  project  in  view  of  published  work  on  performance 
of  wavelet-based  communication  protocols  by  other  researchers. 

Task  3  led  to  introduction  and  evaluation  of  an  algorithm  for  detection  of  cyclostationary  signals. 
This  detector  was  used  to  evaluate  of  the  detectability  of  scale-division  multiple  access  signals  by 
cyclostationary  methods. 
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4  Accomplishments  /  New  Findings 

At  the  inception  of  this  project,  only  a  handful  of  researchers  had  looked  into  the  possibility  of 
using  wavelets  [8,  44,  49,  50]  or  fractal  signals  [57,  58]  in  the  modulation  of  communication  signals. 
The  research  group  at  Arizona  State  University  was  the  first  to  investigate  multiple  access  issues 
in  this  context  [8].  There  are  now  over  forty  published  papers  in  this  area  [2,  5,  9-1.3,  17-21,  24-30, 
33-36,  38-45,  47,  49-50,  53,  55-59]. 

Work  under  this  effort  focused  on  wavelet-based  multiple-access  techniques  in  which  digital 
communication  signals  are  encoded  on  orthogonal  wavelet  sets,  as  was  introduced  in  [8]  and  described 
in  the  original  proposal.  Highlights  of  the  results  obtained  in  this  project  are  summarized  below. 

4.1  Orthonormal  Wavelet  Symbols 

Scale-division  multiple  access  (SDMA)  refers  to  a  multiple-access  communication  scheme  in 
which  the  orthogonal  symbols  are  obtained  by  dilation  and  time  shifting  of  a  single  prototype 
wavelet.  The  various  users’  messages  are  separated  into  channels  based  on  the  scale  of  dilation  of 
the  wavelet  symbol  on  which  their  messages  are  encoded  -  hence  the  term  SDM.4  [8,  10.  52].  .4s 
shown  in  these  references,  the  spectral  structure  of  a  transmitted  SDMA  signal  depends  on  the 
choice  of  the  prototype  wavelet  symbol  in  a  straightforward  way. 

The  goals  in  construction  of  wavelets  for  use  as  symbols  in  SDM.4  are  different  from  those  en¬ 
countered  in  designing  wavelets  for  other  applications.  For  use  in  time-frequency  analysis  of  signals, 
for  example,  it  is  generally  desirable  for  wavelets  to  have  their  energy  simultaneously  concentrated 
in  both  the  time  and  frequency  domains  in  order  to  provide  signal  analysis  that  is  localized  in  the 
phase  plane.  In  spread  spectrum  applications,  however,  wavelets  should  be  compactly  supported  or 
highly  concentrated  in  the  time  domain  while  having  broadband  frequency  structure. 

In  general,  SDMA  relies  on  exact  synchronization  of  aU  transmitters  in  the  network  to  maintain 
orthogonality  of  the  signals  and  avoid  channel  crosstalk.  Even  if  this  is  achieved,  time  delays 
introduced  by  propagation  through  the  transmission  medium  will  cause  the  received  signals  to  be 
non-orthogonal.  This  issue  may  be  addressed  by  requiring  that  the  channels  in  the  network  use 
symbols  with  Fourier  transforms  having  disjoint  supports.  In  this  project,  orthonormal  wavelet 
bases  generated  by  mother  wavelets  whose  frequency  spectra  are  supported  in  several  disjoint  bands 
were  constructed. 

The  simplest  example  of  a  dyadic  orthonormal  basis  of  bandlimited  wavelets  is  generated  by  the 
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mother  wavelet  defined  in  the  frequency  domain  by 


h(u;) 


1  TT  <  |u|  <  27r 
0  otherwise 


(1) 


with  Co  =  2  and  6o  =  1-  A  related,  but  more  complicated,  example  of  an  orthonormal  basis  of 
bandlimited  wavelets  was  described  by  Mallat  using  the  mother  wavelet  defined  by 


h{uj)  = 


1  |a;|  €  (|7r,7r)U[47r,^7r) 

0  otherwise 


(2) 


Examination  of  the  relationship  between  these  two  wavelet  bases  gives  insight  about  how  to  construct 
other  bandlimited  wavelet  bases  in  which  the  dilations  of  the  mother  wavelet  h  do  not  overlap  in 
the  frequency  domain. 

In  what  follows,  the  connected  components  of  the  support  of  a  mother  wavelet’s  Fourier  transform 
h  will  be  called  “slots.”  Throughout  this  section  h  will  be  a  real-valued  wavelet,  which  implies  that 
h(-w)  =  h*(u;)  and  hence  the  slots  are  symmetric  about  w  =  0.  Accordingly,  when  the  values  or 
support  of  h  are  specified  on  the  positive  frequency  a.xis,  they  will  be  assumed  to  be  given  on  the 
negative  frequency  axis  by  symmetry.  The  class  of  wavelets  to  be  constructed  will  be  categorized 
according  to  the  number  of  slots  they  have  in  the  positive  frequency  axis.  The  wavelet  defined  by 
(1)  is  thus  a  “one-slot”  wavelet  while  the  one  defined  by  (2)  is  a  “two-slot”  wavelet.  The  following 
sections  describe  the  construction  of  “n-slot”  orthonormal  wavelet  bases. 

Suppose  the  support  S'o  of  h  is  the  union  of  n  disjoint  intervals 


So  =  [.So,  Si]  U  [S2,  S3]  U  •  •  •  U  [Sn-l ,  Sn]  (.3) 

Then  the  support  of  is 

Sm  =  [.so2~"^,5T2-"^]  U  [s22-"^,.S32-”^]  U  •  •  •  U  [s„_i2-”‘,s„2-"*] 

If  the  measure  of  5^  n5TO'  is  zero  for  all  non-equal  integers  m  and  m'.  So  wiU  be  called  an  orthogonal 
support.  In  particular,  an  orthogonal  support  consisting  of  n  disjoint  intervals  will  be  called  an  n-slot 
orthogonal  support. 

If  ||h||  =  1  and  the  support  of  h  is  orthogonal, 

=  ^(c.A;i.)  =  o 

for  aU  m  7^  m'.  In  this  case,  for  the  time-shifted  and  dilated  replicates  of  h  to  form  an  orthonor¬ 
mal  set  in  £2,  it  remains  only  to  ensure  that  /i"  and  fi™',  are  orthonormal  for  n  n'  and  all  m. 
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This  is  simplified  by  the  observation  that 

{hi,hi)  =  {hi,h^) 

Thus,  one  must  only  verify  that  all  of  the  time-shifted  replicates  of  h  are  orthonormal  at  a  single 
level  of  dilation. 

4.1.1  Construction  of  Orthogonal  Supports 
Consider  rfeal  numbers  0  <  cq  <  •  •  •  <  c„  =  2co.  Then 

[co,  2co]  =  [co,  Cl]  U  [ci ,  C2]  U  [c2,  C3]  U  •  •  •  U  [c„_i ,  2co] 

Let  pi, . . .  be  distinct  integers  and  dilate  each  sub-interval  [cj,  cy+i]  by  2^j+i  to  form  a  set 

5o  =  [co2P‘ ,  ci2?’> ]  U  [ci2P2 ,  C22P2]  u  •  •  •  U  [c„_i2P" ,  2co2P-] 

=  ['50,5l]U  [52,33]  U  •••U  [S2n,52n-l] 

Then  Sq  is  clearly  an  n-slot  orthogonal  support.  Conversely,  suppose  So  is  an  7i-slot  orthogonal 
support  of  the  form  (3).  Then  it  is  straightforward  to  show  that  Sq  can  be  generated  by  the  method 
just  described  with  possibly  some  subintervals  deleted. 

4.1.2  Orthonormality  of  the  Wavelets 

Suppose  ||h||  =  1  and  the  support  of  h  is  orthogonal.  Then  {h^y  is  an  orthonormal  wavelet  set 
if  and  only  if  {hg,  /jq  )  =  0  for  all  integers  n  and  n'  with  n  ^  n'.  This  is  equivalent  to 

0  =  {ho,hQ)=  f  h{t  -  n)h‘{t  —  n')  dt 

for  aU  n  n' .  With  k  —  n  —  n'  ^  this  expression  becomes 

/  1^(0;) 1 2  do;  =  0  (4) 

for  any  non- zero  integer  k.  The  following  section  presents  some  orthonormal  wavelet  sets  corre¬ 
sponding  to  certain  solutions  of  this  equation. 
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4.1.3  Flat  n-Slot  Orthonormal  Wavelets 

The  first  n-slot  orthonormal  wavelets  constructed  were  similar  to  the  wavelets  defined  by  (1) 
and  (2)  in  that  their  Fourier  transforms  are  constant  on  the  n  slots  and  zero  elsewhere.  These  will 
be  called  flat  n-slot  wavelets.  Construction  of  "non-flat”  n-slot  orthonormal  wavelets  is  undertaken 
in  the  next  section  of  this  report. 

Let 

V  =  [co2P‘ , Ci2Pi]  U  [ci2P2, C22P^]  U  •  •  •  U  [c„_i2P",  c„2P'>] 
where  0  <  Cq  <  . . .  <  Cn  =  2co  and  P\,P2i  •  ■  -  ipn  are  integers  and  define 


a 


i=i 


-1/2 


Then  if  h(a;)  —  a  for  u  ^  V  and  h[uj)  —  0  otherwise,  ||/i||  =  1  and  equation  (4)  reduces  to 


i 

y-coZPi 

=  a\  -P 

yci2Pl 

r-ciP2 

/  + 

/  +  /  +• 

-'co2Pi  J 

J  C\2P2 

yci2Pl 

/■C22P2 

rCn2P^ 

=  2o2(  /  + 

/  ’ 

••+  /  )cos(A;; 

Jco2Pi 

/ci2P2 

JCn-\2Pn 

2ofl 

=  -j^{- sin(co2^^Ar) -1- sin 

{ci2^^k)  -  sin(ci2^- 

sin(c„_i2P"/t)  -t-  sin(c„2P"A:)}  =  0 


for  aU  A:  7^  0.  This  is  equivalent  to 


-  sin(co2P‘A;)  -1-  sin(ci2PiA;)  -  sin(ci2PU-)  +  sm(c22^^k) - sin(c„_i2P''A;)  -|-  sin(c„2P"A:)  =  0  (5) 

for  aU  Ar  7^  0.  If  co2pi,  Ci2p^,  Ci2p^,  C22P2,  . . .,  c„_i2P",  and  c„2P"  satisfy  equation  (.5),  then  so  wiU 
A/co2Pi,  Mci2Pi,  iV/ci2P2,  Mc22p^,  . . .,  A/c„_i2P’',  and  Mc„2P"  for  any  natural  number  M.  Thus, 
there  exists  a  support  of  lowest  dilation  level  which  will  be  called  the  mother  support.  In  other 
words,  if  the  end  points  of  the  support 


V  =  [co2P' , Cx2Pi]  U  [ci2P% C22P2]  u  •  •  •  U  [cn_i2P",  c„2P'*] 
satisfy  (5)  and  the  end  points  of  the  dilated  support 

[co2PVM, ci2P7M]  U  [ci2PVM, C22^yM]  U  •  •  -U  [c„_i2P"/M, c„2P7M] 
do  not  satisfy  (5)  for  any  M  G  N,  then  V  is  the  mother  support. 


Example  1  In  the  one- slot  case,  the  support  consists  of  only  one  interval  [co'2^^  ,'2co2^^]  with  cq  >  0. 
Denoting  a  =  co2p*,  equation  (5)  becomes  —  sin(aA:)  -|-  sin(2afc)  =  0,  or 


sin(a&)[2  cos(ak)  —  1]  =  0 

The  only  non-trivial  solution  of  this  equation  that  holds  for  any  non-zero  integer  ifc  is  a  =  m-  with 
m  €  N.  Hence,  the  support  is  [m7r,2m7r],  the  mother  support  is  [7r,27r],  and  an  orthonormal  wavelet 
basis  obtained  having  the  mother  wavelet  defined  in  equation  (1). 

Example  2  In  the  two-slot  case,  the  cut  points  cq  and  ci  and  the  dilation  p  should  satisfy 

-  sin(coA:)  -|-  sin(ciA;)  -  sin(ci2^A;)  -f-  sin(2co2PA;)  =  0 

This  equation  can  be  decomposed  into  two  equations  (with  loss  of  some  solutions) 

{sin(2co2PA;)  -  sin(co^)  =  0 
sin(ci2*’A;)  —  sin(ciAr)  =  0 

which  can  be  solved  simultaneously  under  the  additional  conditions  ||/i||  =  1  and  a  =  1  to  yield 

{2P 

^0  —  2P+ 1  - 1  ^ 

Cl  =  TT 


With  p  =  0  and  cq  =  t,  the  wavelet  obtained  is  the  one  defined  in  equation  (1);  p  =  2  and  cq  =  =- 
yields  the  wavelet  of  equation  (2).  For  p  =  1  and  cq  =  |x,  a  new  bandUmited  orthonormal  mother 
wavelet  defined  by 


1  |u;|  6  [|7r,7r)U  [27r,|7r) 

0  otherwise 


is  obtained. 


4.1.4  Non-Flat  n-Slot  Orthonormal  Wavelets 
Denoting  ^  equation  (4)  becomes 

f  e^'^^\h{uj)\^du=  I  doj  =  2Tnl^{k)  =  0. 

for  aU  non-zero  integers  k.  This  formulation  leads  to  the  following: 


(6) 


Theorem  1  Let  ■0’  be  a  real-valued  and  non-negative  function  with  orthogonal  support.  Further 
suppose  du)  =  27r  and  i){k)  =  0  for  all  non-zero  integers  k.  Then  with  |fi|  =  and  the 

phase  of  h  chosen  arbitrarily  ( except  for  conjugate  symmetry),  h  is  an  orthonormal  mother  wavelet. 
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In  view  of  this  result,  construction  of  the  desired  non-flat  wavelets  hinges  on  finding  functions 
'll)  having  the  properties  specified.  Such  functions  may  be  constructed  as  follows. 

Denote  sinc(^)  =  — and  observe  that  sinc(A:)  =  0  for  all  non-zero  integers  k.  Its  Fourier 
transform  has  values 

^  J  1  lu|  <  X 

sinc(a;)  =  < 

0  otherwise 

Suppose  S  —  [co2Pi,ci2P']  U  •••U  [cn-i2^’' ,  is  an  n-slot  orthogonal  support  in  which  each 

slot  has  measure  larger  than  2x.  Let  /  be  any  non-trivial  function  supported  in  [cq2'^'-  +x,  Ci2^^  - 
x]  U  •  •  •  U  [Cn-i2*’"  +  X,  Cn2P"  —  x]  (figure  1)  and  whose  convolution  with  sine  is  non-negative.  Then 
V’(0  =  /(O  sinc(t)  is  zero  for  all  non-zero  integer  values  of  t,  =  [sine  +  /]  is  non-negative,  and  the 
support  of  ^  is  orthogonal.  Hence,  an  orthonormal  mother  wavelet  h  can  be  obtained  by  letting 


and  setting  the  phase  of  h  arbitrarily. 


\h(' 


2TTi> 


4  du 


(a) 


c  2 
i-/ 


(b) 


c  2 
A'-/ 


c  2"*-;r 
k 


(u 


Figure  1:  (a)  Each  interval  in  the  original  n-band  orthogonal  support  S  has  length  larger  than  2x. 
(b)  The  support  of  /  is  formed  by  shrinking  each  interval  in  5  by  x  on  its  left  and  right  ends. 


Note  that  h  may  be  a  dilated  version  of  another  orthonormal  mother  wavelet.  If  after  rescaling, 
there  exists  some  N  €  N,  such  that  h{2~^t)  satisfies  equation  (4),  i.e.. 


for  aU  nonzero  integers  k,  and 


[  €^^'^\h{2^u)\'^  doj  =  0 

Jl. 


for  some  nonzero  integer  k,  then  h'  defined  by  h'{t)  =  2  ^l‘^h{2  ^t)  is  also  an  orthonormal  mother 
wavelet. 


7 


The  first  part  of  the  following  result  is  a  consequence  of  the  above  construction;  the  converse 
portion  is  proven  in  [10]. 

Theorem  2  For  any  n-slot  orthogonal  support  So,  there  exist  orthonormal  mother  wavelets  having 
support  Sq.  Conversely,  any  n-slot  orthonormal  mother  wavelet  can  be  constructed  in  the  way 
discussed  above. 


Example  3  Let  cq  =  7r,pi  =  1,  ci  =  Ttt/L,  and  C2  =  2.  Then  the  support  of  h  is  [T,77r/4]u[77r,87r]. 
Expand  this  support  by  a  factor  of  4  to  [47r,77r]  U  [287r,32T]  to  make  the  length  of  each  interval 
larger  than  27r.  Choose  the  support  of  /  to  be  [otCjOtt]  U  [297r,.31x]  and  generate  /  subject  to  the 
condition  that  /  *  sine  must  be  real  and  nonnegative;  here  define  /  by 


Then 


1  brr  <  |a;|  <  Ott  and  297r  <  |u;|  <  SItt 
0  elsewhere 


=  [f*  sinc](i^’) 


and  is  as  shown  in  figure  2.  Taking  the  square  root  of  tp  as  the  absolute  value  of  h  and  imposing 
zero  phase,  the  orthonormal  mother  wavelet  with  Fourier  transform  depicted  in  figure  3  is  obtained. 


M/(w) 


Figure  2:  Convolution  of  /  with  the  Fourier  transform  of  a  sine  function  ensures  'tp{k)  =  0  for  all 
non-zero  integers  k. 


4.2  Orthonormal  Wavelet  Symbols  with  Frequency  Overlap 

If  bandlimited  wavelets  are  to  be  used  as  orthogonal  symbols  in  a  low-probability  of  exploitation 
(LPE)  communication  system,  it  may  be  desirable  to  have  the  bands  overlap  to  avoid  the  possibility 
of  an  unintended  receiver  separating  the  channels  using  a  relatively  simple  bank  of  bandpass  filters. 
Following  a  technique  developed  by  Suter  and  Oxley  [46],  bandlimited  wavelet  symbols  have  been 
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Figure  3:  An  orthogonal  mother  wavelet  h  is  obtained  as  the  inverse  Fourier  transform  of  the  square 
root  of  'll;. 


developed  which  (i)  generate  orthonormal  sets  and  (ii)  whose  replicates  at  different  scales  overlap 
in  the  frequency  domain.  The  approach  is  summarized  below. 

The  construction  given  in  [46]  assumes  that  the  real  line  has  been  partitioned  into  disjoint 
intervals  Ij,  j  G  S.  If  {fj,k\k  E  Z}  is  an  orthonormal  basis  for  then  €  Z}  is  an 

orthonormal  basis  for  X^(E),  A  new  orthonormal  basis  is  constructed  by  extension  of  the  functions 
fj^k  off  the  interval  /j  in  a  specific  way,  resulting  in  a  basis  in  which  the  elements  ‘'overlap. 

This  construction  is  modified  to  the  wavelet  case  to  construct  “overlapped”  bandUmited  or¬ 
thonormal  wavelet  symbols  as  follows.  Let  h  be  a  bandlimited  wavelet  having  the  properties: 

1.  The  Fourier  transform  h  of  h  is  real- valued  and  even.  The  steps  in  the  construction  described 
below  apply  to  the  portion  of  h  on  the  positive  frequency  axis  and  are  assumed  to  also  be 
applied  to  the  portion  of  h  on  the  negative  frequency  axis  to  maintain  even  symmetry. 

2.  The  support  of  h  is  orthogonal  in  the  sense  introduced  above;  i.e.,  the  portion  on  the  positive 
frequency  axis  is  of  the  form 

U  •  ‘-U  [2^"an_i,2^-27r] 

where  n  >  1,  tt  <  cti  <  •  •  •  <  <  27r,  and  dj,  ...,dn  are  integers.  The  construction  given 

here  also  assumes  di  ^  2dn^ 

For  integers  j  and  k,  define 

hj,k{t)  =  2-2/1  ^ j  (7) 

and  denote  by  hj^k  the  Fourier  transform  of  hj^k- 

Suppose  {hj^k\j,  k  £  Z,}  is  an  orthonormal  basis  of  X^(M)  (some  examples  of  such  bases  are 
constructed  above  and  others  are  described  in  [6]).  With  e  <  min{(ai  -  7r)/2,(27r  -  a„_i)/2}, 
define  an  extension  h  of  h  by 
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1.  Constructing  the  odd  extension  of  h  about  the  point  2^^7r  into  the  interval  [2^^(7r  —  6),2^^7r]; 

2.  Constructing  the  even  extension  of  A  about  the  point  2^”27r  into  the  interval  [2^^27r,  2^”(27r+6)]; 
and 

3.  Repeating  these  steps  on  the  negative  frequency  axis  to  preserve  even  symmetry. 

An  example  is  shown  in  figure  4  using  Mallat’s  wavelet  defined  in  equation  (1). 

Denote  the  extended  wavelet  just  constructed  by  h  and  let  w  denote  an  frequency  domain 
window  function  with  support  identical  to  that  of  h  and  with  the  amplitude  normalization  properties 
described  in  [46].  An  example  of  w  is  depicted  in  figure  4.  With  u  defined  by 

fi(u;)  =  w{uj)h((jj)  (8) 

the  set  {uj^klJik  G  E}  is  shown  to  be  an  orthonormal  wavelet  basis  of  Also,  by  construe* 

tion,  the  dilated  replicates  of  u  have  non-trivial  overlap.  Figure  5  shows  a  particular  overlapped 
bandlimited  orthonormal  w'avelet  symbol  u{u)  constructed  in  the  above  example. 


•so  5  10  -10  -5  0  S  10 


Figure  4:  Initial  frequency  extension  h(oj)  of  the  wavelet  h  (left)  and  the  window  function  iZ-(a;) 
applied  to  it  (right). 

4.3  Countermeasures 

A  detector  for  cyclostationary  signals  was  developed  during  this  elfort.  This  detector  (depicted 
in  figure  6)  uses  magnitude-squared  coherence  (MSC)  estimation  as  a  measure  of  the  spectral  cor¬ 
relation  in  a  signal.  The  use  of  MSC  estimation  allows  the  statistical  behavior  of  the  detector 
in  a  noise-only  environment  to  determined  analytically  and  threshold  values  corresponding  to  de¬ 
sired  false  alarm  probabilities  to  be  computed  [1-5,  16].  These  capabilities  represent  a  substantial 
improvement  over  similar  cyclostationary  detectors  described  in  the  research  hterature. 
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Figure  5:  Example  of  an  overlapped  bandlimited  orthonormal  wavelet  symbol  u{uj). 


Figure  6:  Discrete-time  version  of  cyclostationary  feature  detector. 

Receiver  operating  characteristic  (ROC)  curves  for  this  detector  against  various  types  of  cy¬ 
clostationary  communication  signals  have  been  determined  by  (i)  using  the  analytical  results  to 
set  detection  a  threshold  for  a  particular  false  alarm  probability,  and  then  (ii)  running  computer 
simulations  with  signal  present  to  estimate  the  corresponding  detection  probability.  Performance 
results  were  obtained  by  applying  the  detector  against  standard  communication  signals  (e.g.,  polar 
BPSK)  and  against  some  LPE-type  communication  signals,  including  direct  sequence  spread  spec¬ 
trum,  transmitted  reference  spread  spectrum,  and  SDMA  signals.  The  case  of  SDMA  signals  was 
of  particular  interest  in  this  project,  and  is  discussed  further  below. 

Consider  an  M-channel  SDMA  signal  r  described  by 

A/— 1  oo 

’■(')  =  E  E  (9) 

771=0  A:=— OO 

where  bm,k  represents  the  bit  on  the  channel  and  is  time-shifted  replicate  of  the 
wavelet  symbol  if;  at  dilation  level  m.  The  spectral  correlation  density  (SCD)  of  such  a  SDMA  signal 
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was  found  to  be 


s:if)  = 


0 


a 


n2^ 


JL 

To 

-iO 


(10) 


where  ip  is  the  Fourier  transform  of  /  is  the  spectral  frequency,  a  is  the  cyclic  frequency,  and  To 
is  the  fundamental  bit  period  (i.e.,  the  time  shift  required  for  orthonormality  of  the  wavelet  symbols 
at  the  lowest  scale  of  dilation).  This  expression  shows  that  the  SCD  of  a  SDMA  signals  will  be 
non-zero  only  at  cyclic  frequencies  that  are  integer  multiples  of  1/To.  This  behavior,  which  has  also 
been  verified  empirically,  indicates  that  SDMA  signals  exhibit  cyclostationary  characteristics  that 
can  possibly  be  exploited  for  the  purpose  of  signal  detection. 

The  cyclostationary  feature  detector  was  applied  to  detect  the  presence  of  a  SDMA  signal  in 
noise.  A  six-channel  SDMA  signal  was  used  as  the  input  signal- for  the  single-cycle  detector.  The 


wavelet  symbol  used  in  the  SDMA  signal  was  the  Daubechies-4  wavelet.  An  initial  analysis  was 
performed  to  determine  at  which  values  of  a  the  spectral  correlation  density  would  be  nonzero.  It 
was  found  that  the  SCD  had  a  strong  component  at  a  =  2^/To  and  /  =  0.  The  single-cycle  detector 


was  set  to  operate  at  these  values  of  /  and  a.  The  number  N  of  independent  samples  used  in  the 
correlations  performed  by  the  single-cycle  detector  was  set  to  32.  The  bandwidth  of  the  single-sided 
bandpass  filters  was  set  to  1.28/To.  Simulations  were  run  at  three  different  SNR’s:  OdB.  3dB  and 
9dB.  One  thousand  outputs  w^ere  used  in  obtaining  the  empirical  cumulative  distribution  for  each 
SNR  setting.  The  ROC  curves  obtained  by  this  procedure  are  shown  in  figure  7.  These  results 
suggest  that  SDMA  signals  are  vulnerable,  at  least  to  some  extent,  to  detection  by  cyclostationary 
methods. 

Another  detector  for  polycyclic  signals  based  on  generalized  coherence  estimation  [7]  has  been 
defined  and  its  application  against  various  types  of  communications  waveforms,  including  SDMA, 
is  currently  under  investigation. 


4.4  Generalized  Wavelet  Transforms 

Work  in  connection  with  this  project  has  led  to  consideration  of  how  the  continuous  wavelet 
transform  can  be  generalized  to  higher  dimensions.  Some  mathematical  results  arising  from  this 
work  and  a  connection  with  image  analysis  are  summarized  in  this  section. 

Several  multidimensional  generalizations  of  the  one-dimensional  wavelet  transform  are  currently. 
In  the  most  widely  known  (and  obvious)  generalization,  one-dimensional  transforms  are  applied 
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Figure  7;  ROC  curves  for  the  single-cycle  detector  with  a  six-channel  SDM.4  signal  at  its  input. 

separately  in  each  orthogonal  direction.  In  what  follows,  this  will  be  referred  to  as  the  rectilinear 
wavelet  transform;  it  is  the  method  used  by  most  wavelet-based  image  compression  techniques. 
Another  wavelet  generalization  is  the  two  dimensional  transform  by  Murenzl[37],  literature  on  which 
can  be  more  easily  found  in  the  paper  by  Antoine  et  al.  [1].  Here,  it  will  be  called  the  circular 
wavelet  transform,  since  it  generalizes  the  dilations  of  the  wavelet  transform  to  the  set  of  dilations 
and  rotations  on  Neither  of  these  generalizations  encompasses  the  other.  Thus,  a  question  of 
interest  is  whether  there  is  a  stronger  generalization  that  encompasses  both.  The  remainder  of  this 
section  describes  a  generalized  wavelet  transform  that  includes  as  special  cases  the  one- dimensional, 
rectilinear,  and  circular  wavelet  transforms.  .4mong  the  multidimensional  wavelet  transforms  that 
arise  from  this  approach  are  several  that  appear  to  be  unknown  in  engineering  applications.  The 
circular  wavelet  transform  is  shown  to  arise  as  a  result  of  the  relationship  between  and  C  in  the 
context  of  this  generalization  and  is  further  shown  to  be  strongly  related  to  the  cortex  transform 
of  Watson  [54].  An  implementation  of  the  cortex  transform  from  the  circular  wavelet  transform  is 
developed  as  an  example  of  this  relationship. 

4.4.1  Theoretical  development 

Recall  first  the  essentials  of  one-dimensional  wavelet  analysis,  details  of  which  can  be  found  in 
[14].  With  a,  5,  X  €  E,  A  €  X^(E),  and  h  the  Fourier  transform  of  h,  define 


1.3 


and 


Ch 


Ij 


|A(x)|2 


dx. 


ya 

If  Ch  is  finite,  h  is  admissible.  In  this  case,  the  reconstruction  theorem  states  that  any  /  G  X^(E) 
can  be  reconstructed  from  its  wavelet  coefficients  ^/,  by  the  formula 


The  rectilinear  wavelet  transform  generalizes  the  one-dimensional  transform  to  a  separable  mul¬ 
tidimensional  transform  by  letting  a,6,x  G  IR",  i.e., 


a  = 

ctl 

b  = 

■  b^  ' 

X  = 

Xi 

with  the  following  alterations  in  the  definitions  and  theorem: 

(x)  = 


\ 

1 

ai 

^ _ h 

1 

\JCl\  •  •  •  dn 

n 

/ 

Ch=  [  dx 

-’Xn\ 

Ck  7lR^7ln  aj  .  \  / 

This  generalization  is  equivalent  to  taking  the  one-dimensional  wavelet  transform  successively  in  each 
orthogonal  direction  of  the  multidimensional  space.  For  n  =  1,  the  rectilinear  wavelet  transform 
reduces  to  the  one-dimensional  wavelet  transform. 

An  alternative  generalization  of  the  one-dimensional  wavelet  transform  to  is  the  circular 
wavelet  transform  [1,  37],  for  which  a  G  E  is  replaced  with  the  pair  (a,0)  G  E  X  [0,27r)  and  6  G  E  is 
replaced  by  6  =  [61,62]'^  €  E^.  Let  rg  denote  the  matrix  which  performs  rotation  in  E^  by  9  degrees. 
Then  the  definitions  and  theorem  for  circular  wavelet  analysis  are 

c.=  / 

f=-^f  ry  adbdOda. 

Ch  Jr  Jo  Jr^  ^  ' 


14 


Clearly,  the  circular  wavelet  transform  has  a  different  character  than  the  rectilinear  wavelet 
transform  on  In  particular,  it  allows  non-separable  wavelets  in  two  dimensions  and  replaces 
independent  dilation  in  each  dimension  by  rotation  and  dilation  simultaneously  in  both  dimensions. 
Despite  their  distinct  characters,  both  the  rectilinear  and  circular  wavelet  transforms  arise  as  special 
cases  of  a  more  general  multidimensional  wavelet  transform,  the  formulation  of  which  is  based  on 
ring  multiplications  on  E^. 

Consider  E’^  as  a  vector  space  over  the  field  E.  If  a  ring  multiplication  ‘*0”  is  introduced  to  E”' 
so  that  E"^  becomes  an  algebra  under  the  multiplication  then  there  is,  for  each  a  €  E^,  a  left  regular 
representation  :  E’^  ->  E’^  defined  by  LaX  —  a  0  x.  There  is  also  an  associated  linear  operator 
5a:  :  E’^  ^  E^  defined  by  Sxd  =  Suppose  the  ring  multiplication  has  the  property  that  the  set 

of  elements  with  ring  inverses  are  dense  in  E’^.  Then,  with  h  6  X^(E^)  and  a,  6,  x  £  E.^,  the  following 
definitions  and  theorem  form  the  mathematical  foundation  of  a  generalized  multidimensional  wavelet 
analysis: 


(x)  = 


1 


Ch=  [ 


\/|  det  La 

\h(xt- 


det  5^ 


h{La-\x-b)), 
dx. 


As  examples  of  this,  suppose  the  ring  multiplication  on  E”  is  defined  by 


ai 

■  bx  ' 

axbi 

* 

0 

z= 

.  bn  . 

then  La  and  Sa  are  the  diagonal  matrices  defined  by  La  =  Sa  =  diag(ai, a2,  •  •  • ,  «n)-  This  makes  the 
generalized  definitions  and  theorem  equivalent  to  the  definitions  and  theorem  of  rectilinear  wavelet 
analysis.  Since  rectilinear  wavelet  analysis  with  n  =  1  is  the  same  as  one-dimensional  wavelet 
analysis,  both  one-dimensional  wavelet  analysis  and  rectilinear  wavelet  analysis  arise  as  special 
cases  of  this  formulation. 

Similarly,  circular  wavelet  analysis  comes  from  the  ring  multiplication  “o”  on  defined  by 


ai 

0 

'  bx  ' 

(iib\  —  a2^2 

_  0.2 

.  ^2 

0\b2  "h  02b\ 

This  ring  multiplication  is  same  as  the  multiplication  of  the  complex  numbers.  In  this  case 


and  the  generalized  definitions  and  reconstruction  theorem  become  equivalent  to  the  definitions  and 
reconstruction  theorem  of  circular  wavelet  analysis. 


4.4.2  Connection  to  the  cortex  transform 

Circular  wavelet  analysis  has  fundamental  variables  of  scale  (which  corresponds  to  logarithmic 
frequency),  rotation,  and  two-dimensional  positional  shift.  These  correspond  to  the  basic  organiza¬ 
tion  of  Watson  s  cortex  transform[54],  which  was  devised  as  a  means  to  process  images  in  a  similar 
way  to  the  processing  in  the  human  visual  cortex.  Hence  there  is  a  relationship  between  the  circular 
wavelet  transform,  the  cortex  transform,  and  image  processing  in  the  human  visual  cortex.  It  is 
shown  here  that  this  relationship  can  be  strengthened  by  actually  deriving  the  cortex  transform 
from  the  circular  wavelet  transform. 

The  cortex  transform  was  implemented  by  Watson  by  partitioning  the  frequency  domain  into 
separate  angularly  oriented  regions  at  logarithmically  spaced  frequency  steps,  bandlimiting  the 
image  to  each  region,  and  sampling  the  bandlimited  pieces  as  efficiently  as  possible  using  available 
techniques  from  sampling  theory. 

Following  a  similar  procedure,  the  frequency  domain  can  be  divided  according  to  the  figure 
below,  where  zero  frequency  is  in  the  center: 


This  division  corresponds  to  a  discrete  sampling  of  the  frequency  and  rotation  variables  of  the 
circular  wavelet  transform.  In  order  to  bandlimit  the  image  to  a  given  region,  the  ‘‘mother  wavelet” 
h  is  chosen  so  that  its  Fourier  transform  h  is  the  characteristic  function  a  bin,  which  consists  of  two 
opposing  frequency  regions  as  shown  below: 
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Once  the  image  is  bandlimited  to  a  bin,  positional  samples  are  used  to  represent  the  bandlimited 
function.  These  positional  samples  produce  frequency  replications.  Aliasing  is  avoided  by  preventing 
these  frequency  replications  from  overlapping.  An  example  of  an  optimal  sampling  for  a  bin  is  shown 
below. 


This  frequency  replication  pattern  is  according  to  theory  the  optimally  efficient  replication  that 
avoids  aliasing,  but  in  practice  aliasing  is  not  avoided  due  to  filter  imperfections.  Thus,  it  is  often 
desirable  to  introduce  a  controlled  amount  of  inefficiency  in  order  to  reduce  aliasing.  In  this  case, 
the  frequency  replication  pattern  looks  more  like  this: 


The  efficiency  of  this  sampling  is  controlled,  so  it  is  easy  to  make  it  more  efficient  than  the  fixed 
efficiency  sampling  proposed  by  Watson. 

Preliminary  results  show  that  the  transform  adequately  represents  an  image  and  reconstructs 
it  from  its  representation  (figure  8).  In  order  to  obtain  these  results,  special  consideration  must 
be  paid  to  preserving  frequencies  near  zero,  since  the  transform  has  a  singularity  there.  It  is  also 
crucial  that  the  image  be  of  relatively  large  size  (at  least  256  x  256  pixels),  since  the  processing 
assumes  an  image  is  continuous,  and  the  violation  of  this  assumption  by  a  coarsely  sampled  image 
introduces  significant  artifacts. 
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Figure  8:  The  Leena  image  (left)  and  the  result  after  it  has  been  decomposed  and  reconstructed  by 
the  transform  (right). 
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[6b-2]  A,  Clausen  and  D.  Cochran,  “An  Invariance  Property  of  the  Generalized  Coherence  Esti¬ 
mate,”  To  appear  in:  IEEE  Transactions  on  Signal  Processing. 
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6c.  Recent  submissions 

[6c-l]  A.  Clausen  and  D.  Cochran,  ‘‘Generalized  Coherence:  Invariance  of  Estimate  to  Distri- 
bution  of  One  Channel,”  Submitted  to  the  30th  Asilomar  Conference  on  Signals,  Systems,  and 
Computers,  November  1996. 

[6c-2]  L.L.  Scharf  and  D.  Cochran,  “On  the  Distribution  of  Linear  Prediction  Error,”  Submitted 
to  the  30th  Asilomar  Conference  on  Signals,  Systems,  and  Computers,  November  1996. 

[6c-3]  J.Q.  Trelewicz  and  D.  Cochran,  “A  Sampling  Approach  to  Bandlimited  Instantaneously 
Companded  Signals,”  Submitted  to  the  30th  Asilomar  Conference  on  Signals,  Systems,  and  Com¬ 
puters,  November  1996. 

6d.  Other  related  publications 

[6d-l]  D.  Cochran,  “Review  of  Cyclostationarity  in  Communications  and  Signal  Processing," 
(Invited  review)  in  IEEE  Communications  Magazine,  vol.  32,  no.  4,  April  1994.  Also  in  Proceedings 
of  the  IEEE,  October  1994. 

7  Interactions/Transitions 

7a.  Meetings,  Conferences,  and  Seminars 

Numerous  university  and  industrial  colloquium  presentations  have  been  presented  on  the  research 
performed  in  this  project.  These  include: 

•  D.  Cochran  presented  a  talk  entitled  "Wavelet -Based  Multiple- Access  Spread  Spectrum”  in  the 
Motorola  Government  Systems  and  Technologies  Group  Signal  Processing  Series  [16  attendees] 
in  March  1993. 

•  D.  Cochran  presented  a  talk  entitled  “Detection  of  Cyclostationary  and  Polycyclostationary 
Signals  at  the  DSP  group  meeting  at  MIT  [17  attendees]  in  October  1994.  During  the  same 
visit,  he  met  with  A.  Willsky,  C.  Karl,  and  R.  Learned  to  exchange  ideas  on  wavelet  and 
wavelet-packet  methods  in  spread  spectrum. 

•  D.  Cochran  presented  a  talk  entitled  “Multidimensional  Wavelet  Generalizations  and  Watson’s 
Cortex  Transform”  in  the  Vision  Seminar  at  Harvard  [25  attendees]  in  September  1994. 

•  D.  Cochran  presented  a  talk  entitled  “Wavelet  Applications  in  Communications”  in  a  seminar 
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sponsored  jointly  by  the  DSP  and  Communications  groups  at  Georgia  Tech  [45  attendees]  in 
October  1994. 

In  addition,  D.  Cochran  participated  in  the  following  conference  activities  related  to  the  topic 
of  this  project: 

•  He  organized  and  chaired  an  invited  session  entitled  “Applications  of  Wavelets  and  Chaotic 
Signals  in  Communications”  at  the  Asilomar  Conference  on  Signals,  Systems,  and  Comput¬ 
ers  in  November  1993,  Participants  included  R.  Orr  (Atlantic  Aerospace),  S.  IsabeUe  (MIT), 
M.  Motamed  and  A.  Zakhor  (UC  Berkeley),  C.  Wei  (Arizona  State  University),  and  J.S.  Gold¬ 
stein  (USAF  and  Georgia  Tech). 

•  He  w'as  a  keynote  speaker  (with  J.J.  Benedetto,  M,V.  Wickerhauser,  and  W.  Sweldens)  at  a 
wavelet  applications  workshop  in  Melbourne,  Australia  in  February  1995.  The  topic  of  this 
presentation  was  w’avelet  applications  in  communication  systems. 

•  He  is  scheduled  to  be  plenary  speaker  (with  L.  Cohen)  at  the  ANZIIS-96  Conference  in  Novem¬ 
ber  1996. 

7b.  Consultive  and  Advisory  Functions 

During  this  project,  D.  Cochran  consulted  informally  with  J.  Stephens  at  Wright  Laboratory 
and  with  Capt.  J.S.  Goldstein  and  Dr,  A.  Lindsey  of  Rome  Laboratory  on  applications  of  wavelets 
and  cyclostationary  signal  processing  in  covert  communications  and  countermeasures.  He  has  also 
been  collaborating  with  the  Centre  for  Signal  and  Information  Processing  Research  (CSSIP)  and 
the  Defence  Science  and  Technology  Organization  (DSTO)  in  Australia  on  wavelet  applications  in 
communications  and  other  topics. 

7c.  Transitions 

None. 


8  New  Discoveries 

No  inventions  or  patent  disclosures  have  resulted  from  this  research  program. 
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9  Honors/ Awards 

None. 
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